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^ ■ ABSTRACT 

I'T ^ ! We have recently discovered a supermassive binary black hole system with 

ly-j ! a projected separ ation between the two black holes of 7.3 parsecs in the radio 



o 

u 



^ ! galaxy 0402+379 ( Rodriguez et al.ll2006l ). This is the most compact supermassive 

binary black hole pair yet imaged by more than two orders of magnitude. We 
present Global VLBI observations at 1.3464 GHz of this radio galaxy, taken to 
improve the quality of the HI data. Two absorption lines are found toward the 



P_i! southern jet of the source, one redshifted by 370 ± 10 km s^^ and the other 

Q i blueshifted by 700 ± 10 km s'^ with respect to the systemic velocity of the 

^ i source, which, along with the results obtained for the opacity distribution over 

ci '. the source, suggests the presence of two mass clumps rotating around the central 
region of the source. We propose a model consisting of a geometrically thick disk, 

^ , of which we only see a couple of clumps, that reproduces the velocities measured 

I from the HI absorption profiles. These clumps rotate in circular Keplerian orbits 

^ I around an axis that crosses one of the supermassive black holes of the binary 

I system in 0402+379. We find an upper limit for the inclination angle of the twin 

^ I jets of the source to the line of sight of 6* = 66°, which, according to the proposed 

Q\ \ model, implies a lower limit on the central mass of ~ 7 x 10® Mq and a lower 

^ ' limit for the scale height of the thick disk of ~ 12 pc . 
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1. Introduction 

Given tliat most galaxies harbor supermassive black holes at their centers (Richstone 
et al. 1998; Gebhardt et al. 2000), and that galaxy mergers are common, binary black 
holes should likewise be common. An understanding of the evolution and formation of these 
systems is important for an understanding of the evolution and formation of galaxies in 
general (Silk & Rees 1998, Merritt 2006). Theoretical descriptions of supermassive binary 
black hole systems and their accretion disks have been investigated by Hayasaki et al. 2007, 
2008 and by MacFadyen & Milosavljevic 2008. 

Our ability to resolve the supermassive black holes in any given binary system depends 
on the separation between them, on their distance from Earth, and on the resolving power of 
the telescope used. It is believed that the longest timescale in the evolution of a supermassive 
binary black hole system leading up to coalescence is the stage in which the system is closely 
bound (~ 0.1 - 10 pc), meaning that in most of these systems the black hole pair can only 
be resolved by VLBI observations (in the case where both black holes are radio loud) which 
provides resolutions of milliarcseconds and finer. This c ould explain why very few such 



systems have been found (see review by iKomossal (j2003al ) detailing observational evidence 
for supermassive black hole binaries). 

Some source properties like X-shaped radio galaxies and double-double radio galax- 
ies, helical radio-jets, double-horned emission line profiles, and semi-periodic variations in 
lightcurves have been taken as indirect evidence for compact binary black holes though other 
explanations are possible. Some wider systems have, however, been found more directly. For 
example, the ultra luminous galaxy NGC 6240, discovered by the Chandra X-ray observa- 



tory, w as found to have a pair of active supermassive black holes at its center (IKomossa et al. 



2003bl ). separated by a distance of 1.4 kpc. Another system that has been known for some 



time is the double AGN (7 kpc separation) constitutin g the radio source 3C 75, which was 



discovered by the VLA to have two pairs of radio jets (jOwen et al.lll985l ). 



The radio galaxy 0402+379 was recently found to contain two central, compact, flat 
spectrum, variable components (designated CI and C2), with a projected separation of 7.3 
pc, a feature which had not been observed in any other compact source, making this system 
the most compact supermassive binary black hole pair yet imaged by more than two orders of 
magnitude, with an estimated system mass of a few 10^ Mq (Maness et al. 2004; Rodriguez 
et al. 2006). 



Maness et al.l (120041 ) performed spectral line observations of 0402+379 at 1.348 GHz. 
HI absorption at a redshift of 560 km s~^ from the systemic velocity of the source (16,489 
± 300 km s~^; Xu et al. 1994) was observed and attributed to a high velocity gas system. 
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possibly due to a merger. Following t his discovery, a spectr um of the HI region was taken in 



2004 with the Westerbork telescope (IMorganti et al 



20091 ) . which showed two components 



with velocities separated by 1000 km s~^. In this article we present Global VLBI observa- 
tions taken to improve the quality of the HI data on 0402+379, with the purpose of better 
understanding the origin of the deviation from the systemic velocity seen in the HI gas in 
this source. 



At the redshift of 0402+379 of 0.055, Ho=75 km s^^ Mpc'^ and qo 
mas = 1.06 pc is obtained. 



0.5, a scale of 1 



2. Observations 

2.1. Global VLBI Observations 

Global VLBI observation^ were made on 2007 March 17 at 1.3464 GHz. A single 
intermediate frequency with a bandwidth of 16 MHz was observed with 256 channels in 
both right and left circular polarizations, resulting in a frequency resolution of 62.5 kHz, 
corresponding to a velocity resolution of 15 km s~^. Four level quantization was employed. 
The net integration time on 0402+379 was 497 minutes. 

Standard flagging, amplitude calibration, fringe-fitting, and bandpass calibration (3C 
111 was used for both gain and bandpass calibration) were followed in the Astronomical 
I mage Processing System (AIPS; van Moorsel et al. 1996). AIPS reduction scripts described 



m 



Ulvestad et al.l (120011 ) were used for a large part of the reduction. Spectral line Doppler 



corrections were also applied in AIPS. All manual editing, frequency averaging procedures 



imaging, deconvolution, and self-calibration were done using Difmap (jShepherd et al.lll995l ). 
A clean cube was also produced by first performing a continuum subtraction using the task 
UVLSF in AIPS, then using Difmap in order to obtain a cleaned map for each spectral 
frequency, and finally combining all the maps into a cube with the task MCUBE in AIPS. 



^Global VLBI — Very Long Baseline Array + Effelsberg + Jodrell Bank + Westerbork + Green Bank 
Telescope + Onsala. 
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Results 



3.1. Radio Continuum 



Figure [T] shows a naturally weighted 1.3 GHz image of 0402+379 from the 2007 Global 
VLBI observations. The source consists of two diametrically opposed jets and a central 
region containing two active nuclei not resolved at this frequency (see Rodriguez et al. 2006 
for a detailed multi- frequency study of this radio galaxy). The northern jet is pointing in 
the northeast direction whereas the southern jet is pointing in the southwest direction. The 
orientation of t he source at this fre quency is consistent with that se en by the VLA at bot h 
1.5 and 5 GHz jManess et al.ll2004l ) and by the VLBA at 0.3 GHz jRodriguez et al.ll2006[ l. 
The image shown in panel (a) was tapered and restored with a circular 25 mas synthesized 
beam, in order to show better the extended structure in this radio galaxy. The source spans 
~ 1000 mas (~ 1000 pc). The image shown in panel (b) was tapered and restored with a 
8.15 X 3.72 mas synthesized beam. We see structure on scales of ~ 500 mas (~ 500 pc). 



3.2. HI Absorption 

Figure [2] shows HI absorption profiles taken from four regions of the source. The contin- 
uum has been subtracted from the spectra using the task UVLSF in AIPS, which removes a 
continuum model from the u,v data of all channels. The contours are taken from the 2007 
Global VLBI observations at 1.3 GHz and the color scale from the 2005 VLBA observations 
at 5 GHz. Two absorption lines are evident, which appear to be at two different locations 
toward the southern jet. From this point on we will designate CW the western component, 
where we see the stronger line; and CE the eastern component, where we see the weaker line. 
Measurements of these lines are given in Table [H where we also show the peak opacity and 
the column density for both components, which we calculated according to, 

/oo 
T{v)dv(km s"^) 
-oo 

~ 1.8224 X 10^%(K)Ai;(km s'^) ^ t{v), 

(1) 

where Ts is the spin temperature, Av is the velocity resolution, and the summation is over 
the channels where we detect absorption for CW and CE respectively. 

We determined the central velocity of the HI absorption lines to be 16,927 ± 7 km s~^ 
and 15,856 ± 9 km s^^ for CW and CE respectively. Fr om the most recent rep orted redshift 



for 0402+379, as measured from optical emission lines (jRodriguez et al.l 120061 ). we find that 
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the observed line for CW is redshifted by 370 ± 10 km s~^ and for CE is blueshifted by 700 
± 10 km s~^ from the systemic velocity of the source (16,558 ± 3 km s~^). No absorption 
was found at either CI or C2 (see Figure [2]). We calculate a limit on the HI opacity at the 
location of CI and C2, and obtain < 0.17 and < 0.03 repectively. Figure [3] shows a map of 
both the central velocity and width of the HI absorption profiles over the source. 

Figure m shows a velocity slice of the continuum-subtracted cube (right panel), accom- 
panied by the HI opacity distribution over the source (left panel). This result shows that 
either the two locations where we find absorption are localized regions, rather than being 
part of a more extended and perhaps homogeneous structure, or our sensitivity prevented 
us from detecting HI absorption over a broader region. In order to explore this question we 
calculated the lower limit on the peak opacity across the source assuming an intensity for the 
line of 3(7, where a is the rms noise in a single channel. The result is shown in Figure O We 
see that outside the region where we detect absorption (compare with Figure H]) we would 
need a peak opacity of at least ~ 0.02, comparable to the peak opacity found for components 
CE and CW, in order to have a 3(T detection. The fact that we did not detect this means 
that the peak opacity in this region is even smaller than the lower limit found, suggesting 
that what we are observing are two different, localized clumps, rather than a more extended 
and homogeneous structure. The fact that we do not see absorption against C2 supports 
this scenario. 

We can estimate the HI mass of the clumps from the calculated column densities (Table 
[T]) and assuming a spherical shape, according to the following relation. 



m 



—r NnimH 



(2) 



where itlh is the hydrogen atom mass and r is the radius of the clump estimated from the 
opacity distribution map shown in Figure HI For CW we used r ~ 12 pc and for CE r ~ 8 
pc. We find a mass ranging from ~ (7 — 400) x 10^ Mq for component CW, and ~ (1 — 
60) X IO^Mq for component CE, assuming spin temperatures ranging from 100 K to 6000 
K. These m asses are bigger than typical values of HI cloud masses of 60 Mq found in the 
Milky Way Jstil et al.1120061 ). 



4. Discussion 

The results found from our 2007 Global VLBI observations at 1.3 GHz show two HI 
absorption lines, one blueshifted and the other redshifted with respect to the systemic ve- 
locity. The redshifted line (component CW) shows a FWHM of 300 ± 20 km s~^, and the 
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blueshifted line (component CE) a FWHM of 170 ± 20 km s~^. This absorption must be 
from neutral hydrogen along the line-of-sight to the central continuum source, but this gas 
could be close to the nucleus or far away; could be falling in, outflowing, or in rotation. 

Assuming a virialized cloud of ~ 10^ Mq and a radius of 5 pc we get, from the virial 
theorem, a velocity for the H atoms of ~ 3 km s~^. On the other hand, if we use the high 
end temperature of 10^ K the thermal speed of an H atom is ~ 15 km s~^. Thus, under these 
assumptions, at most ~ 20 km s~^ of the line widths observed for both CE and CW could 
be accounted for, leaving a considerable balance which could be attributed to the proximity 
of these components to the central parsecs of the source, where the gravitational potential 
well is deepest. 

Another possibility is that the broad line widths and large velo cities for one or both 



components are due to a jet-cloud interaction. lOosterloo et al.l (l2000l ) studied the Seyfert 2 
galaxy IC 5063, where a strong interaction between the radio jet and a molecular cloud of the 
interstellar medium (ISM) is ocurring at the position of the western radio lobe of the source. 
The most prominent absorption feature found is blueshifted over 600 km s~^ with respect 
to the systemic velocity, a value that lies outside the range allowed by rotational kinematics 
of the large-scale HI disk. As an effect of the outflow produced by the interaction as well 
as the geometry of the source, there will be components moving both away and towards the 
observer. Since only absorption in those components in front of the radio continuum can be 
observed, the blueshifted component was the only one detected. Returning to our case, it is 
very unlikely that the motion found in both components, CW (redshifted component) and CE 
(blueshifted component), which are only a few parsecs apart, is due to jet-cloud interactions, 
since this would require a dramatic change in the way the southern jet is moving, on scales 
of only a few parsecs. 

We consider that both CW and CE are part of the same rotating disk structure, of which 
we only see a couple of clumps, and for the remainder of the discussion, we assume that this 
is the case and propose a model that reproduces the observations. We also assume that the 
mass of component CI is significantly smaller than that of component C2, thus not affecting 
the stability of the disk. This type of binary, where the two black holes have different masses, 
are thought to occur more frequently, since minor galactic mergers are more common tha n 



major mergers in hierarchical models of galaxy formation (lArmitage fc Natarajan 1 12002| ). 



We realize that in a thi n disk the binary would likely cle ar out all material within twice 



the radius of the binary (IMacFadyen fc Milosavljevidl2008l ). but perhaps this process is less 



efficient for a thick disk, or material is resupplied. 
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4.1. Rotating Disk 

If we assume that components CW and CE are rotating about C2 in circular Keplerian 
orbits lying on the same plane, the only way of reproducing the measured velocities (accord- 
ing to our proposed model, explained below) is by requiri ng a very high central mass, of at 



least ~ IO^^Mq, in disagreement with previous estimates (IRodriguez et al.ll2006l ) as well as 
typical supermassive black hole masses. 

Thus, we propose a model consisting of a geometrically thick disk, rather than a ge- 
ometrically thin disk, rotating in a Keplerian circular orbit, which can be part of a larger 
scale torus, that can reproduce the velocities that we observe from the HI absorption, at the 
measured locations. We assume that CW and CE are rotating in parallel planes, both offset 
from the plane in which C2 lies. The distances between these planes are free parameters in 
our model. We call dl the separation between C2 and the center of the circular orbit 1, that 
CE is following; and d2 the separation between C2 and the center of the circular orbit 2, 
that CW is following. Figure E] is a cartoon showing the configuration of the model we are 
proposing. 

We assume that t he rotation axis crosses component C2, which is at the systemic velocity 



of 16,558 ± 3 km s~^ (IRodriguez et al.ll2006l ). a natural assumption considering that the jets 
of this source seem to be emerging from C2, not from CI. Also, we suggest that the two 
absorption lines that are observed could be explained by the presence of two mass clumps, 
both part of a clumpy geometrically thick disk structure (an idea that is supported by the 
opacity distribution map obtained. Figure H]). 

We let dl and d2 vary over a range of values and calculate, for each different combination 
{dl,d2) what the orientation in space of the thick disk should be in order to best reproduce 
the observed velocities, via a minimization routine. After obtaining the inclination of the 
thick disk for each pair {dl,d2), we again apply a minimization routine in order to find the 
combination {dl,d2) that best reproduces the measured velocities. 

This same procedure was applied with different values of the central mass, which allowed 
us to set a lower limit on this quantity of ~ 1.4 x 10^ Mq (by requiring that the difference 
between the velocities predicted by the model and the observed velocities was smaller than 
la). See Appendix Rl for more details on the model. There is no upper limit on the mass 
that we can find using this model, since for greater masses we can always find an inclination 
that reproduces very accurately the measured values (the inclination angle, with respect to 
the line of sight, decreases as we increase the central mass). 



Rodriguez et al.l (120061 ) made an estimate of the central mass of 0402+379, based on an 



optical spectrum of the core of this radio galaxy taken with the Hobby-Eberly Telescope, 
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which shows a red shoulder suggesting two components with a velocity separation of ~300 
km s~^. At the observed projected separation between the nuclei of 7.3 pc, an orbital velocity 
of 300 km s~^ implies a system mass of 1.5 x 10^ Mq. For this particular value of the central 
mass, our model gives an inclination of the rotation axis of ~ 85° with respect to the line 
of sight, thus requiring the source to be very close to the plane of the sky. A value for d2 
of ~ 12 pc is obtained, which sets a lower limit on the scale height of the thick disk. This 
value is consistent with observed parameters of supposed tori (see for example Peck et al. 
1999, where a 20 pc scale height torus is estimated for the radio source 1946+708). 

Considering the size of our synthesized beam and the extension of the two regions where 
absorption is observed, we were not able to measure a significant velocity gradient within CW 
or CE. The small variations in velocity found are not reproduced by our model. However, 
other factors could influence these mass concentrations that we are not taken into account 
in our model, such as proper motions associated with each clump. 



4.2. More Constraints on the Inclination of the Disk 

Assuming that C2 is the origin of the radio emission on parsec scales, we can constrain 
the orientation of 0402+379. In the simple beaming model for simultaneously ejected jet 
components moving in opposite directions, the arm length ratio D dep ends on the intrinsic 



speed [5 = v/c and the angle of the twin jets to the line of sight 9 (ITaylor fc Vermeulen 



19971 ). The arm length ratio, -D, is given by 



z^ = ^=fi±^i (3) 



do, \1 — (3 cos 9 

where the apparent projected distances from C2 (assumed to be the origin of radio emission) 
are for the northern jet (a pproaching side) and (is for the southern jet (receding side). 



From this we get (3 cos 9 = 0.4 (IRodriguez et al.ll2006l ). a result that implies that the intrinsic 



velocity must be at least 0.4c and 9 must be less than 66°. Based on the observed morphology 
of the source, we can assume a minimum value for ^ of ~ 30°. According to our proposed 
model (see §4.11) these constraints on the inclination set a lower limit on the central mass of 
the system of ~ 7 x 10^ Mq, for 9 = 66°. Angles as small as ~ 30° are allowed, however, 
they require a very large mass (~ 10^° 



Figure [7] is a cartoon showing the results obtained by our model for a central mass of 
7 X 10^ Mq {9 = 66°), overlaid with 5 GHz contours of 0402+379 from the 2005 VLBA 
observations. We can see that the rotation axis aligns very well, within ~ 10°, with the 
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jet axis, which was not a constraint imposed by the model. If we asume that the center of 
the orbit that CI follows is located at the position of C2, and also that component C2 is 
significantly more massive than CI, we obtain an orbital velocity for CI of ~ 580 km s^^ 
(for a mass of 7 x 10^ Mq and 6 = 66°). This is consistent with the current limit from VLBI 
monitoring of < 26, 400 km s~^. 

5. Conclusion 

Global VLBI observations at 1.3 GHz were performed on the radio galaxy 0402+379, 
the most compact supermassive binary black hole pair yet imaged. Two absorption lines 
were found toward the southern jet of the source, one redshifted by 370 ± 10 km s^^ and the 
other blueshifted by 700 ± 10 km s~^ with respect to the systemic velocity of the source. 

A model consisting of a thick disk, of which we only see a couple of clumps, was de- 
voloped in order to reproduce the velocities measured from the HI absorption profiles. These 
clumps, components CW and CE, rotate in circular Keplerian orbits around an axis that 
crosses one of the supermassive black holes of the binary system in 0402+379, component 
C2. We found an upper limit for the inclination angle of the twin jets to the line of sight of 
9 — 66°, which, according to the proposed model, implies a lower limit on the central mass 
of ~ 7 X 10^ Mq and a lower hmit for the scale height of the thick disk of ~ 12 pc . 

Facilities: Very Long Baseline Array, Effelsberg, Jodrell Bank, Westerbork, Green Bank 
Telescope, Onsala 25m telescope 
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A. Rotating Disk Model 

Figure [6] is a cartoon stiowing tlie configuration of tfie model we are proposing. We show 
component C2, which is chosen as the origin, as well as components CE, with coordinates 
(xl,yl), and CW, with coordinates (x2,y2). The two circular orbits drawn are the trajectories 
CE and CW would follow, respectively, on the rotating thick disk. Also, we show the 
projection of the system on the plane of the sky, since that is what we actually observe. 

Let us assume our disk has an arbitrary orientation in space, specified by the normal 
unitary vector n = {nx,ny,nz)- Since n is unitary we have, 

nl + nl + nl = l (Al) 

According to Kepler's Laws, a body moving around a mass M, at a radius R, has a 
velocity given by. 

Thus, the angular velocity is given by, 

00 = VGM{Ry/^ (A3) 

If To is the angular velocity vector of our rotating thick disk, and R is any point on the 
disk, then the tangential velocity V is given by, 

V = 1j xR (A4) 



where uj = oon and R = r — d , where d is the separation between the origin and 
the center of the circular orbit we are considering, and 1^ is the position vector of the point 
we are considering (in our case, the coordinates of either CE or CW). Thus, we have, 

R = {Rx, Ry, Rz) = {x + drix, y + duy, z + driz) (A5) 

Let the plane x-y be the plane of the sky and z be the line of sight direction. We want 
to know the z component of the tangential velocity V , because that component is the one 
we can measure (from our absorption lines), which is given by: 

Vz = V .z = uin^Ry - UyR^) (A6) 

Finally, we can also use the fact that Zj and R have to be orthogonal, that is, 

Ij .R = =^ rixX + Uyy + n^z + d = (A7) 
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We can now put everything together and get an expression for \4, 

Vz = VGM [n^{y + duy) - ny{x + dn^)] [{x + dn^f + {y + duyf + {z + dn^f] (A8) 

where 

^_-d-xn^-yn^ and = J I - nl - nl (A9) 



Thus, we have an equation that gives us the z component of the tangential velocity of 
any point with a circular orbit in our rotating disk, in terms of the x and y components of 
the normal unitary vector of the disk, and Uy, and in terms of the x and y components 
of any point in the disk. 

There are two locations in our source, where we can measure the velocity. Let us define 
Vzi^o and Vz2,o as the observed velocities of CE and CW respectively; and Vzi^c and Vz2,c the 
calculated velocities according to Equation IA8I Finally let us define the function that we will 
minimize in order to obtain the inclination of the thick disk that reproduces the measured 
velocities best, 

T_ l^zl.o ~ Kl,c| \Vz2,o — Vz2,c\ f\-\n\ 

~ \v \ \v \ ^ ^ 

For a particular value of the central mass of the system, M, we went through the 
following routine in order to find the inclination of the thick disk, more precisely and ny 
(since Hz can be determined in terms of these two quantities): We calculated the function 
J for every combination of (n^., Uy), with the constraint that n is unitary. By minimizing 
J we found the best pair {n^, Uy). Since dl and d2 are also free parameters, we let them 
vary over a range, so for every combination {dl, d2) we followed the procedure mentioned 
above and found the best pair (n^, ny). After having done this for every combination {dl, 
d2), again by a minimization routine, we found the best pair {dl, d2). Thus, for a particular 
value of the central mass we were able to find the inclination of the disk that reproduced 
the observed velocities best. This same routine was also followed for different values of the 
central mass, which let us find a lower limit on this quantity of ~ 1.4 x 10^ Mq, by requiring 
that the difference between the velocities predicted by the model and the observed velocities 
was smaller than la, with: 



-2 + ^2 (All) 



{Vzl, o-Vzl,c? {Vz2,o - Vz2,cf 

where Ai and A2 are the errors in the measurements of Vzi^o and Vz2,o respectively. 
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Table 1. Gaussian Functions Fitted to HI Absorption Profiles in the Southern Jet. 



Component * 


Amplitude (mjy) 


Central Velocity (km s ^) 


FWHM (km s'^) 




Nhi (cm 2) 


CW 


2.8 ± 0.1 


16,927 ± 7 


300 ± 20 


0.025 ± 0.001 


(1.303 ± 0.006) xlO^i » 












(7.82 ± 0.03) xl022 b 


CE 


1.5 ± 0.2 


15,856 ± 9 


170 ± 20 


0.018 ± 0.002 


(4.1 ± 0.1) xlO^" 












(2.44 ± 0.07) xl022 b 



* C W and CE refer to the western and eastern jet components where we find absorption lines respectively. 
^Assuming a spin temperature of 100 K 
^Assuming a spin temperature of 6000 K 
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Fig. 1.— Naturally weighted 2007 Global VLBI image of 0402+379 at 1.3 GHz. Contours 
are drawn beginning at 3a and increase by factors of 2 thereafter. In (b) the synthesized 
beam is 8.15 x 3.72 mas (shown in the bottom left corner). The peak flux density is 0.13 
Jy/beam and rms noise is 0.02 mJy. In (a) the image was tapered and restored with a 
circular 25 mas synthesized beam (shown in the bottom left corner). The peak flux density 
is 0.46 Jy/beam and rms noise is 0.05 mJy 
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Fig. 2. — HI absorption profiles taken from four regions of the source (components CI, C2, 
C E, and CW). The arrow in the center of each spectrum shows where the systemic velocity 
is (IRodriguez et al.ll2006h . and the velocity resolution is 15 km s~^. The rms noise in a single 
channel is 0.14 mJy/beam. The contours are taken from the 2007 Global VLBI observations 
at 1.3 GHz and are set at lOcr, increasing by a factor of 2 thereafter; and the color scale 
image is taken from the 2005 VLBA observations at 5 GHz. 
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Fig. 3. — (a) Map of the central velocity and (b) width of the HI absorption profiles in the 
source. Both maps were generated fitting Gaussian functions at each pixel where at least a 
3cr detection of a line was found. The color scale in both images is in units of km s~^. The 
contours are taken from the 2007 Global VLBI observations at 1.3 GHz and are set at 3a, 
increasing by a factor of 2 thereafter. 
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Fig. 4. — Velocity slice of the continuum-subtracted cube described in §3.21 (right panel), ac- 
companied by the integrated HI opacity distribution over the source (left panel) , in which we 
indicate where the slice was taken. The map was generated by combining a continuum im- 
age of the source with the continuum-subtracted cube. To generate the opacity distribution 
shown we did a one-dimensional zeroth moment fitting to each row, requiring a continuum 
emission of at least 0.12 mJy and a line emission of at least 0.75 mJy (corresponding to ~ Sa 
respectively). The contours are taken from the 2007 Global VLBI observations at 1.3 GHz 
and are set at 3a, increasing by a factor of 2 thereafter. 
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Fig. 5. — Calculated lower limit on the peak opacity across the source assuming an intensity 
for the line of So". The contours are taken from the 2007 Global VLBI observations at 1.3 
GHz and are set at Sa, increasing by a factor of 2 thereafter. 
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Fig. 6. — Configuration of tlie model we are proposing. Sliown is component C2, cliosen as 
tlie origin, as well as components CE, witli coordinates (xl,yl), and CW, witli coordinates 
(x2,y2). Tlie two circular orbits drawn are the trajectories CE and CW would follow, respec- 
tively, on the rotating thick disk, shown in yellow. An inclination angle of ~ 75° between 
the rotation axis and the line of sight was used. 
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Fig. 7. — Results obtained from the rotating thick disk model, projected on the x-y plane 
(RA and DEC plane), overlaid with 5 GHz contours from the 2005 VLBA observations. The 
blue section of the rotation axis is pointing towards us, whereas the red section is pointing 
away. For a mass of 7 x 10^ Mq an inclination angle of ~ 66° between the rotation axis and 
the line of sight is obtained. The position of component C2, chosen as the origin, is shown, 
as well as the position of CI, CE and CW. The two green circular orbits drawn are the 
trajectories CE and CW would follow, respectively, on the rotating thick disk. The yellow 
circular orbit is the trajectory CI would follow assuming C2 as the center of the orbit. 



